Vibrational Tamm states at the edges of graphene nanoribbons 
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We study vibrational states localized at the edges of graphene nanoribbons. Such surface oscil- 
lations can be considered as a phonon analog of Tamm states well known in the electronic theory. 
We consider both armchair and zigzag graphene stripes and demonstrate that surface modes corre- 
spond to phonons localized at the edges of the graphene nanoribbon, and they can be classified as 
in-plane and out-of-plane modes. In addition, in armchair nanoribbons anharmonic edge modes can 
experience longitudinal localization in the form of self- localized nonlinear modes, or surface breather 
solitons. 



I. INTRODUCTION 

Graphite nanocrystallitcs are recently discovered 
carbon-based nanomaterials 1 which are a few atoms thick 
and stable under ambient conditions, and they keep many 
promises as materials for nanoscale transistors 2 . Re- 
markable properties of graphite structures make their 
study one of the hot topics of nanoscience 3 . 

Graphene nanoribbons are effectively low-dimensional 
nanostructures similar to carbon nanotubes, but they 
differ from nanotubes by the presence of edges. Due 
to the edges, graphene nanoribbons can demonstrate 
many novel geometry-driven properties, depending on 
their width and helicity. A majority of the current stud- 
ies of graphene nanoribbons target their electronic and 
magnetic properties stipulated by the presence of edges, 
including the existence of the edge modes 4 , which are an 
analog of surface states in this two-dimensional geometry. 

In general, surface states are spatially localized modes 
which can be generated at surfaces 5 , and such surface 
states have been studied in many branches of of physics 
including electrons in crystals 6 ' 7 , surface phonons 8 , sur- 
face polaritons 9 , and optical surface modes in waveg- 
uide arrays 10 ' 11 . In this paper, we analyze the specific 
properties of vibrational phonon modes localized at the 
edges of graphene nanoribbons, the so-called edge phonon 
modes. Such localized surface modes can be considered 
as a phonon analog of Tamm states well known in the 
electronic theory 6 . 

Vibrational edge modes in graphene nanoribbons have 
been discussed recently for the hydrogen-terminated 
graphene nanoribbons 12 . Using the second-generation 
reactive empirical bond-order potential and density- 
functional theory calculations, it was demonstrated that 
two different edge modes can exist in armchair nanorib- 
bons, whereas the zigzag nanoribbons can support only 
one such mode. Importantly, all those edge modes were 
found to correspond to the out-of-plane vibrations of the 
edge hydrogen atoms. 

In this paper, we employ numerical analysis and molec- 
ular dynamics simulations to study localized modes in 
graphene nanoribbons. We demonstrate that such nanos- 
tructures can support much larger number of vibrational 



phonon modes localized at the edges. More specifically, 
we find that the armchair nanoribbon can support 7 edge 
modes, whereas the zigzag nanoribbon can support 12 
edge modes. Only one part of such modes correspond to 
out-of-plane oscillations of the edge atoms, whereas and 
the second part is characterized by in-plane edge oscilla- 
tions. Chemical modification of the edge atoms can lead 
to the generation of larger number of localized surface 
modes. In addition, numerical simulations of anharmonic 
oscillations demonstrate that in armchair nanoribbons 
the vibrational energy can become self-trapped along the 
edge of the stripe with the formation of both in-plane 
and out-of-plane surface breather solitons. 

The paper is organized as follows. In Sec. II we discuss 
our model including the potentials for describing the in- 
teratomic interactions and the edge effects. Section III is 
devoted to a general linear analysis, and it discusses the 
dispersion curves for both armchair and zigzag nanorib- 
bons. In Sec. IV we describe the specific properties of 
the edge modes and corresponding energy localization. 
Some properties of anharmonic oscillations and surface 
breather solitons are discussed in Sec. V, whereas Sec. VI 
concludes the paper. 



II. MODEL OF GRAPHENE NANORIBBONS 

We model the graphene nanoribbon as a planar stripe 
of graphite, with the properties depending on the stripe 
width and chirality. Here we consider two most typical 
geometries of the armchair nanoribbon [sec Fig. 1(a)] 
and zigzag nanoribbon [sec Fig. 1(b)]. 

The structure of the nanoribbon can be presented 
through a longitudinal repetition of the elementary cell 
composed of K atoms, as shown in Figs. l(a,b). For the 
armchair, the number K is multiple of 4, whereas for the 
zigzag nanoribbon the number K is multiple of 2. We use 
the atom numbering shown in Figs. l(a,b). In this case, 
each carbon atom has a two-component index a = (n, k), 
where n = 0, ±1, ±2, ... stands for the number of the ele- 
mentary cells, and k = 1, 2, K stands for the number 
of atoms in the cell. 

Each elementary cell of the armchair nanoribbon has 
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FIG. 1: Schematic view of (a) armchair and (b) zigzag 
nanoribbons with the atom numbering. The edge atoms are 
shown as filled circles. Dotted lines separate the elementary 
cells of the nanoribbon. M is the number of atoms in the 
elementary cell. 



four edge atoms, whereas the zigzag nanoribbon has only 
two edge atoms. In Figures l(a,b), wc show these edge 
atoms as filled circles. In a realistic case, the edge atoms 
are chemically modified and, generally speaking, the fol- 
lowing situations are possible 13,14 : a single hydrogen CH 
(hydrogen-terminated graphene nanoribbon) , CH2 - two 
hydrogens, COH - a hydroxyl group, and CHOH - a hy- 
drogen atom and a hydroxyl group that resulted from the 
water decomposition. 

In this study, we do not take into account the spe- 
cific chemical nature of the surface atoms, but consider a 
change of the effective mass of the edge atom. Therefore, 
in our model of graphene nanoribbons wc take the mass of 
atoms inside the stripe as Mq = 12m p , and for the edge 
atoms we consider a larger mass, such as Mi = 13m p , 
14m p , 29m p , 30m p (where m p = 1.6603 ■ 10~ 27 kg is the 
proton mass). 

To describe nanoribbon oscillations, we write the sys- 
tem Hamiltonian in the form, 
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M(n, k ) (U( n ,fc) 3 U( n ,fc) ) + P( n ,k) 



, (i) 



where M a is the mass of the hydrogen atom with the 
index a = (n, k) (for internal atoms we take M a = Mq, 
whereas for the edge atoms we take a larger mass, M a = 
All > Mq), u a = (x a (t),y a (t), z a (tj) is the radius-vector 
of the carbon atom with the index a at the moment t. 
The term P a describes the interaction of the atom with 
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FIG. 2: Configurations containing up to i-th nearest-neighbor 
interactions for: (a) i = 1, (b) i = 2, (c) i = 3, (d) i = 4, and 
(e) i = 5. 



the index a = (n, k) and its neighboring atoms. The 
potential depends on variations of bond length, bond an- 
gles, and dihedral angles between the planes formed by 
three neighboring carbon atoms, and it can be written in 
the form 

p = E c/i + E c/2 +E c/ 3+E c/ 4 + E c/ ^ (2) 

where Qj, with i = 1,2,3,4,5, are the sets of configura- 
tions including up to nearest-neighbor interactions. Ow- 
ing to a large redundancy, the sets only need to contain 
configurations of the atoms shown in Fig. 2, including 
their rotated and mirrored versions. 

The potential U\ (u Q , Ug) describes the deformation en- 
ergy due to a direct interaction between pairs of atoms 
with the indexes a and (3, as shown in Fig. 2(a). 
The potential C/2(u a ,u^,u 7 ) describes the deformation 
energy of the angle between the valcnt bonds u Q u^ 
and ugii 7 , see Fig. 2(b). Potentials Ui(u a , u^, u 7 , us), 
i = 3, 4, 5, describes the deformation energy associated 
with a change of the effective angle between the planes 
u Q , Ug, u 7 and u^, u 7 , 115, as shown in Figs. 2(c,d,e). 

We use the potentials employed in the modeling of the 
dynamics of large polymer macromolecules 15-19 for the 
valent bond coupling, 

C/i(ui,u 2 ) = ei{exp[-aoGo-/9o)] - I} 2 , p = |u 2 -ui|, 

(3) 

where t\ = 4.9632 eV is the energy of the valent bond 
and po = 1.418 A is the equilibrium length of the bond; 
the potential of the valent angle 

C/ 2 (ui,u 2 ,u 3 ) = e 2 (cos(^ - cosc^o) 2 , (4) 
costp = (u 3 - u 2 ,ui - u 2 )/(|u 3 - u 2 | • |u 2 - Ui|), 

so that the equilibrium value of the angle is defined as 
costpo = cos(2-7r/3) = —1/2; the potential of the torsion 
angle 

E7i(ui, 112,113,114) = £i(l - z 4 cos0), (5) 
COS0 = (vi, v 2 )/(|vi| • |v 2 |), 
Vl = (U2 - Ul) x (113 - u 2 ), 
v 2 = (u 3 - u 2 ) x (u 3 - u 4 ), 

where the sign z,; = 1 for the indices i = 3, 4 (equilibrium 
value of the torsional angle </>o = 0) and Zi = — 1 for the 
index i = 5 (<fo = 7r). 

The specific values of the parameters are «o = 
1.7889 A" 1 , e 2 = 1.3143 eV, and e 3 = 0.499 eV, and 
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they are found from the frequency spectrum of small- 
amplitude oscillations of a sheet of graphite 20 . Accord- 
ing to the results of Rcf. 21 the energy £4 is close to the 
energy e 3 , whereas £5 -C £4 ( | / £4 1 < 1/20). Therefore, 
in what follows we use the values £4 = £3 = 0.499 eV and 
assume £5=0, the latter means that we omit the last 
term in the sum (2). 



III. DISPERSION CURVES 

In the equilibrium state, {u° n k \}^^l^ k=1 , all atoms 
are in the plane of the nanoribbon, so that all valent 
bonds have the equilibrium length p = po, and all 
angle between the bonds coincide, if — 2n/3. Arm- 
chair nanoribbon is characterized by the longitudinal 
shift h x = Spo and width D y = \/2>(K/2 — l)po/2, 
where K is the number of atoms in the elementary cell. 
Zigzag nanoribbon is characterized by the longitudinal 
shift h x = \/?>pq and width D y = po(3K — 4)/4. 

Next, we introduce 3AT-dimcnsional vector, 

Xn = ( u (n,l) ~ U (n,l)J — ! U {n,K) ~ U °n,if)) 

describing a shift of the atom of the n-th cell from its 
equilibrium position. Then, the nanoribbon Hamiltonian 
can be written in the following form, 




H = ^{-(Mx n ,x„) +'P(x Tl _x,x n ,x n+ i)}, 



(0) 



FIG. 3: (color online) Structure of 3K dispersion curves for 
(a) armchair nanoribbon with the width D y = 11.05 A (K = 
20), and (b) zigzag nanoribbon with the width D y — 
11.34 A (K = 12). Blue (black) curves correspond to the 
in-plane vibrations, and red (gray) curves correspond to out- 
of-plane vibrations. Mass of the edge atoms is Mi = 13m p . 



where M is the diagonal matrix of masses of all atoms of 
the elementary cell. 

Hamiltonian (6) generates the following set of the equa- 
tions of motion: 

-Mx n = 7 , Xl (x„,x n+ i,x n+2 ) -fP X2 (x n _i,x„,x n+ i) 
+7 , X3 (x, l -2,x„_i,x„) (7) 

In the linear approximation, this system takes the form 

-Mx„ = Bix„+B 2 x„ + i+B2X„_i+B 3 x„ + 2+B3X„_2, 

(8) 

where the matrix elements arc defined as 



Bi =P XlXl +V-. 



x 2 x 2 + fx 3 x 3 , 



B2 — "Pxix 2 + "Px 2 X3, B 3 = "Pxix 3 , 

and the matrix of the partial derivatives takes the form 
d 2 V 

-(0,0,0), i,j = 1,2,3. 



1 3 dxidi j 

Solutions of the system of linear equations (8) can be 
sought in the standard form, 



x„ = Aecxp(iqn — icut). 



(9) 



where A is the mode amplitude, e is the normalized di- 
mensionless vector ((Me,e)/Mo = 1), q £ [0,n] is the 



dimcnsionlcss wave number, and u> is the phonon fre- 
quency. Substituting the expression (9) into the system 
(8), we obtain the eigenvalue problem 



w 2 Me = [Bi+B2e* 



B*e 



B 



3e 



2iq_ 



V, 



-2«n 



(10) 



Therefore, in order to find the dispersion relations 
characterizing the modes of the nanoribbon for each fixed 
value of the wave number < q < tt we need to find 
the eigenvalues of the Hcrmitian matrix (10) of the or- 
der 3K x 3K. As a result, the dispersion curves are 
composed of 3K branches {wj(g)}^ 1; as shown in Fig. 
3. Two thirds of the branches corresponds to the atom 
vibrations in the plane of the nanoribbon xy (in-plane 
vibrations), whereas only one third corresponds to the 
vibrations orthogonal to the plane (out-of-plane vibra- 
tions), when the atoms are shifted along the axes z. 

As is seen in Fig. 3, the spectrum of the nanorib- 
bon oscillations occupies the frequency interval [0,w m ], 
where the maximum (cutoff) frequency u> m practically 
does not depend on the chirality. The maximum fre- 
quency u> m = 1600 cm" 1 . This value agrees well with 
the experimental data for a planar graphite 22 ' 23 . Four 
branches of the dispersion curves start from the ori- 
gin, (q = 0, uj = 0). The first branch cui(q) corre- 
sponds to the orthogonal (out-of-plane) bending vibra- 
tions of the nanoribbon; the second branch u>2{q) de- 
scribes the bending planar (in-plane) vibrations in the 
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plane. These branches approach smoothly the axes q: 
UJ (<l)/ ( l — > 0, when q — > 0, so that the corresponding long- 
wave phonons possess no dispersion. The third w 3 (q) and 
fourth W4 (q) branches of the dispersion curves correspond 
to the out-of-plane twisting and in-plane longitudinal vi- 
brational modes, respectively. The corresponding long- 
wave modes possess nonzero dispersion, so that we can 
define two limiting values, 

t ,. w 3 (q) u 4 (q) 
Vt = n x lim , and V\ — h x lim 

q^-Q q g-s-0 q 

which defines the sound speeds of twisting and longitu- 
dinal phonons in the nanoribbon. The values of these 
speeds depend on the nanoribbon chirality. Nanorib- 
bons of the same width but different chiralities support 
phonons with different sound speeds: armchair nanorib- 
bon with K = 20 (width D y = 11.05 A) has the veloci- 
ties Vt = 3400 m/s and u; = 24900 m/s, whereas for the 
zigzag nanoribbon with K = 12 (width D y = 11.34 A) 
these velocities are smaller, v t = 1900 m/s and vi = 
13200 m/s. When the width of the nanoribbon grows, 
the velocity of the longitudinal sound remains unchanged 
while the velocity of the twisting modes decreases mono- 
tonically. 

IV. EDGE MODES 

Solution of the eigenvalue problem (10) demonstrates 
that the nanoribbon can support vibrational modes lo- 
calized at the edges, the so-called edge phonons. To find 
such surface modes, we should consider the nanoribbons 
of a sufficiently large width, and we take the large number 
of atoms in the elementary cell, K = 120. 

To analyze oscillations (9), wc define the distribution 
function of the oscillatory energy in the elementary cell 
of the nanoribbon as follows 

Pl = A/ ? (|e3 ( . 1 _ 1)+1 | 2 + |e 3(l -i) + 2| 2 + |e3( l -i) + 3| 2 )/Mo > 0, 

where i = 1, 2, K is the number of atoms in the cell, 
Mi is the mass of the i-th atoms, ej is the component of 
the eigenvector e = {e^}?^. The energy distribution is 
normalized by the following condition, 

K 

Y,P* = (Me,e)/M = 1. 

i=l 

We introduce the parameter characterizing the transverse 
energy localization in the nanoribbon as follows, 

i=l 

This parameter characterizes the inverse width of the en- 
ergy localization in the elementary cell. If the vibrational 
mode is localized only on one atom (i.e. there exists an 



atom io for which pi = 1), the inverse width is d = 1. 
In the opposite limit, when the vibrational energy is dis- 
tributed equally on all atoms (p, = l/K), we have d = K, 
so that in a general case 1 < d < K. 

To be more specific, we define the vibrational mode as 
being localized provided d < 20. Our numerical analysis 
of the oscillatory cigenmodes shows that localization may 
occur only at the edges of the nanoribbon. The structure 
of the 3K dispersion curves of the armchair and zigzag 
nanoribbons are shown in Figs. 4 and 5, where thick lines 
mark the edge modes of the nanoribbon. Such localized 
surface modes exist in both armchair and zigzag nanorib- 
bons. 

Now we consider in more details the edge modes with 
the wave numbers q = and q = tt. The values of the 
frequencies of these modes are summarized in Table I. 
When the mass of the edge atom is Mi = 13m p (the 
effective edge atom corresponds to the molecular group 
CH), the armchair nanoribbon has 7 edge modes (3 out- 
of-plane and 4 in-planc), and zigzag nanoribbon has 12 
edge modes (6 out-of-plane and 6 in-plane) . If the mass 
of the edge atoms is increased, the number of localized 
modes growths as well. For example, for Mi = 30m p 
(when the effective edge atom corresponds to the molec- 
ular group CHOH), the armchair has already 14 edge 
modes (4 out-of-plane and 10 in-plane), and the zigzag 
nanoribbon has 13 edge modes (6 out-of-plane and 7 in- 
plane) . 

We notice that for K = 120, an increase of the mass of 
the edge atoms does not change substantially the disper- 
sion curves because the number of the edge atoms is much 
smaller than the number of atoms in the elementary cell 
(elementary cell of the armchair nanoribbon has 4 edge 
atoms, and the zigzag nanoribbon - 2 edge atoms). How- 
ever, the surface modes are localized at the edge atoms, 
so that a change of the mass of the edge atoms change 
the mode frequency substantially. As can be seen from 
Figs. 4 and 5, a change of the mass of the edge atoms 
from Mi = 13m p to Mi = 30m p leads mainly to a shift of 
the dispersion curves corresponding to the edge modes. 
This behavior provides an additional indication of the 
existence of strongly localized surface modes with the 
properties different from those of the bulk modes. 

Transverse distribution of the energy of the surface lo- 
calized states in the graphene nanoribbon is shown in 
Figs. 6 and 7. As follows from these results, the oscil- 
lations are mainly localized near the nanoribbon edges. 
For the armchair nanoribbon, the oscillations penetrate 
inside the stripe for less than 4 nm, and for the zigzag 
nanoribbon the oscillations decay for no more than 2 
nm. For Mi = 13m p , the edge atoms of the armchair 
nanoribbon always participate in the surface vibration, 
whereas in zigzag nanoribbons the edge atoms may not 
contribute into the oscillatory mode, as shown in Fig. 7 
for the modes with the numbers I = 6, 7, 8 9). 

The specific structure of the edge oscillations in the 
graphene nanoribbon are shown in Figs. 8 and 9. These 
edge modes can be divided into three groups. The first 
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TABLE I: Frequencies ui (dimension [uj] =cm" 1 ) of the localized edge oscillations with the wave numbers q — and q = it, 
depending on the nanoribbon chirality and mass of the edge atom Mi (dimension [Mi] = m p ). Out-of-plane vibrations are 
marked by italic. Oscillations with the wave number q = are underlined. Corresponding energy densities of the edge modes 
are presented in Figs. 6 and 7. 
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227 1452 102 
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236 268 
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912 
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1478 
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547 1276 172 


284 


368 413 


432 


661 


716 


742 


1136 


1159 


1455 






FIG. 4: (color online) Structure of 3K dispersion curves 
of the armchair nanoribbons for (a) mass of the edge atoms 
Mi = 13m p and (b) Mi = 30m p . The number of atoms in the 
elementary cell is K = 120. Blue (black) curves correspond 
to in-plane vibrations, red (gray) curves - to the out-of-plane 
vibrations, thick black curves - oscillations localized at the 
nanoribbon edges. The marked points show the band edges 
where nonlinear breathers can appear in the form of the lon- 
gitudinally localized edge modes. 



group is composed by low-frequency out-of-plane vibra- 
tions which mainly occur due to deformation of the soft 
dihedral angles. For the armchair nanoribbon with the 
mass of the edge atoms Mi = 13m p , these oscillations 
have the frequencies u> = 102, 150 and 227 cm -1 , whereas 
for the zigzag nanoribbons the similar modes have the 
frequencies: ui = 261, 431, 432, 583, 661, and 716 cm -1 . 
As can be seen from Figs. 6(a), 7(a), and Table I, these 
oscillations are the most localized. The main vibrational 
energy is concentrated in the edge atoms, which oscillate 



FIG. 5: (color online) Structure of 3K dispersion curves of 
the zigzag nanoribbons for (a) mass of the edge atoms Mi = 
13m p and (b) Mi = 30m p . The number of atoms in the 
elementary cell is K = 120. Blue (black) curves correspond 
to in-plane vibrations, red (gray) curves - to the out-of-plane 
vibrations, thick black curves - oscillations localized at the 
nanoribbon edges. 



orthogonally to the nanoribbon plane, except the modes 
in the zigzag nanoribbon with the frequencies 431, 432, 
661, and 716 cm -1 (here the edge atoms are static while 
the hexagons which contains these atoms rotate). 

The second group of the edge modes is composed by 
the edge in-plane vibrations which occur due to a defor- 
mation of the valent angles. For the armchair nanorib- 
bon these modes have the frequencies 236 and 268 cm" 1 , 
and for the zigzag nanoribbon, the frequencies 393 and 
406 cm" 1 . 

Finally, the third group of the surface modes is com- 
posed by high-frequency edge oscillations which occur in 
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(a) 




FIG. 6: Transverse distribution of energy of the edge modes 
in the armchair nanoribbon with the wave numbers q = 0, 
and q = n for the mass of the edge atoms: (a) Mi = 13m p 
(q = for I < 2, q = tt for I > 3) and (b) Mi = 30m p (q = 
for i < 4, g = 7r for i > 5). Here y is the distance to the edge 
of the nanoribbon, I is the mode number, p(y) is the energy 
density normalized by the condition: J p(y)dy = 1. The 
values of the frequencies are summarized in Table I. 



the nanoribbon plane which occur due to the deformation 
of the hard valcnt bonds. For the armchair nanoribbons, 
these are the edge modes with the frequencies 1067 and 
1452 cm -1 , and for the zigzag nanoribbon, with the fre- 
quencies 912, 1194, 1466, and 1478 cm" 1 . 

We would like to mention that in the paper 12 the 
edge vibrations were studied by another method by us- 
ing the second-generation reactive empirical bond order 
(REBOII) potential and density-functional theory calcu- 
lations. This method does not allow to study nanorib- 
bons of large width so that the authors of Ref. 12 consid- 
ered the nanoribbons with K = 20 atoms of hydrogen in 
the elementary cell, taking into account only the oscilla- 
tions of the C H groups. They found several types of the 
out-of-plane edge vibrational modes: two edge modes for 
the armchair nanoribbon (with the frequencies 365 and 
858 cm -1 ), and one mode for the zigzag nanoribbon (with 
the frequency 610 cm -1 ). 

In our work, using the mechanical model and direct 
molecular-dynamics numerical simulations allow to study 
much wider nanoribbons and, therefore, to find many 
more types of surface modes. We consider the edge 
C-H atomic group as a single edge atom with the mass 
Mi = 13m p , so that we are not able to study the vibra- 



(a) 




FIG. 7: Transverse distribution of energy of the edge modes 
in the zigzag nanoribbon with the wave numbers q = 0, and 
q = n for the mass of the edge atoms: (a) Mi = 13m p and (b) 
Mi = 30m p (q — for I < 2, q = n for / > 3). Here y is the 
distance to the edge of the nanoribbon, / is the mode num- 
ber, p(y) is the energy density normalized by the condition: 
Jo P(v)dy ~ 1- The values of the frequencies are summarized 
in Table I. 

tional mode associated with the motion of one hydrogen 
atom with the frequency 858 cm -1 from the paper 12 . 
However, for other two vibrational modes with the fre- 
quencies 365 and 610 cm -1 , both carbon and hydrogen 
atoms participate in vibrations, so that we can find sim- 
ilar modes by applying our approach, with the frequen- 
cies 227 and 583 cm -1 . Taking into account the motion 
of the hydrogen atoms leads to a correction of the os- 
cillation frequency, as well as to a larger number of the 
edge modes, due to additional modes induced by sepa- 
rate motion of hydrogen and carbon atoms but not the 
whole group C-H. 



V. LONGITUDINAL LOCALIZATION AND 
SURFACE SOLITONS 

The surface oscillations described above are localized 
at the edges of the stripe, and they correspond to quasi- 
one-dimensional excitations of the nanoribbon. There- 
fore, similar to other types of nonlinear modes, we may 
expect that such surface mode may experience the longi- 
tudinal localization with the formation of breather soli- 
tons. Such novel type of surface breathers exist due to 
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FIG. 8: Six edge modes of the armchair nanoribbon. Shown are the displacements of atoms of the nanoribbon at quarter of 
an oscillation period, for: (a) out-of-plane vibration with lo = 227 cm - , q — 0; (b) in-plane vibration with u> = 1452 cm -1 , 
q = 0; (c) out-of-plane vibration with ui = 102 cm -1 , q = tt; (d) out-of-plane vibration with uj — 150 cm -1 , q — ty; (e) in-plane 
vibration with co = 236 cm -1 , q = ty; (f) in-plane vibration with uj = 268 cm" 1 , q — ty. Mass of the edge atom is Mi = 13m p . 



the transverse energy localization in the surface state and 
nonlinearity-induced longitudinal localization along the 
nanoribbon edge. 

To analyze the existence of breather solitons, we con- 
sider the nanoribbon of a fixed length and free edges with 
N = 500 periods. We introduce friction at the far edges 
in order to absorb phonons propagating away from the 
localization region. To find spatially localized nonlinear 
states, we integrate numerically the system of equations 
(7) (1 < n < N) with the initial conditions corresponding 



to a localized state at the edge of the nanoribbon, 

x n (i) = A n e cxjpiiqri — iuit), (11) 

where A n = A scch{/i[(n — N/2) + l/2}}, with the param- 
eter /x > characterizing the inverse width of the initial 
excitation, and the amplitude A > 0. A complex unit 
3.A /-dimensional vector e = e r + ie;, a solution of the 
linear equation (10), corresponds to the edge mode with 
the wave numbers q = 0, n and the frequency u>. 

For the system of equations (7), the ansatz (11) corre- 
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FIG. 9: Six edge modes of the zigzag nanoribbon. Shown are the displacements of atoms of the nanoribbon at quarter of an 
oscillation period, for: (a) out-of-plane vibration with u = 583 cm - , q = 0; (b) in-plane vibration with uj — 1466 cm" , q = 0; 
(c) out-of-plane vibration with u = 261 cm -1 , q — ir; (d) in-plane vibration with ui = 393 cm -1 , q = 7r; (e) in-plane vibration 
with uj — 406 cm -1 , q = n; (f) out-of-plane vibration with oj = 431 cm -1 , q = %. Mass of the edge atoms is Mi — 13m p . 



sponds to the following initial conditions, 

x„(0) = A n [e r cos(qh) - eism(qn)}, (12) 
x„(0) = A n uj[e r sm(qn) — cos(gn)], 

where n = n- N/2 + 1/2. 



The breather frequency can only split-off from one of 
the linear dispersion bands. As follows from Figs. 4 and 
5, the frequency branches corresponding to the edge vi- 
brations have either minimal or maximum values at q = 
and q = tt, so that we use these values in the initial con- 
ditions (12). 
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By changing the value of the initial amplitude A and 
the inverse width /x, and analyzing the dynamics of the 
system (7) with the initial conditions (12), we may ar- 
rive at the conclusion if the system under considera- 
tion can support spatially localized states in the form 
of the breather solitons, as self-trapped nonlinear modes 
splitting-off from the corresponding edge of the linear 
dispersion band. If such a breather exists, the initial 
excitation (12) will experience no substantial spreading, 
whereas it will spread out if the breather conditions are 
not satisfied. 

Our extensive numerical simulations reveal that the 
edge breather are not possible for the zigzag nanorib- 
bons, and the initial excitation of the form (12) always 
spreads. However, for the armchair nanoribbons the lon- 
gitudinal localization of the edge modes becomes possi- 
ble with the wave number q = n and frequency u> = 102 
(out-of-plane vibration) and 268 cm -1 (in-plane vibra- 
tion), see Fig. 10. For the inverse width /i = 0.05 and 
amplitude A = 0.008 A, the initial localized vibrational 
state with the frequency uj = 102 cm -1 decays slowly, 
but this decay stops when the initial amplitude is in- 
creased to the value A = 0.016 A, and it evolves into 
a localized breathing mode with the breather frequency 
103 cm -1 . When we increase the initial amplitude to 
the value A = 0.032 A. we observe the generation of two 
breathers which interact nonlinearly. 

Similarly, for the planar edge vibration with the fre- 
quency uj = 268 cm -1 the initial localized state spreads 
when its amplitude is A = 0.004 A [see Fig. 10(a)], it gets 
self-trapped for A = 0.008 A and generates a breather 
with the frequency 267 cm -1 [sec Fig. 10(b)], and for 
A = 0.016 A it creates two interacting breathers. 

Our studies demonstrate that armchair nanoribbons 
can support two types of spatially localized surface 
breathers. One type of breathers originates from the lo- 
calization of out-of-plane edge mode with the frequency 
uj = 102 cm . This edge mode is shown in Fig. 8(c), 
where the edge atoms of the neighboring elementary cells 
oscillate with the opposite phases in the direction per- 
pendicular to the nanoribbon plane. The other type 
of breathers originates from the longitudinal localiza- 
tion of the in-plane edge mode with the frequency uj = 
268 cm -1 . This edge mode is shown in Fig. 8(f), where 
the edge atoms oscillate transversally in the nanoribbon 
plane. The frequency of the former breather grows with 
the amplitude, while the frequency of the latter breather 
decays with the amplitude, in accord with the structure 
of the dispersion curves presented in Fig. 4. 



VI. CONCLUSIONS 

We have demonstrated that both armchair and zigzag 
graphene nanoribbons can support various types of vi- 
brational Tamm states in the form of phonons localized 
at the edges of the graphene stripe. We have shown that 
such localized modes appear in the gaps of the frequency 
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FIG. 10: Dynamics of the initially localized edge excitation 
of the armchair nanoribbon (wave number q = 7r, frequency 
uj = 268 cm -1 , and inverse dimensionless width /j, — 0.05) for 
the initial amplitude: (a) A = 0.004 and (b) A = 0.008 A. 
Shown is the temporal evolution of the the kinetic energy E n 
in two first edge layers of atoms of the nanoribbon (k = 1, 2, 
3, 4). Mass of the edge atoms is Mi = 13m p . 



spectrum of the nanoribbon oscillations, and they can 
be classified as in-plane and out-of-plane surface modes. 
In the former case, all atoms oscillate in the plane of the 
nanoribbon, whereas in the latter case the atoms oscillate 
perpendicular to the plane. We have found that some of 
these edge modes can experience self-localization along 
the nanoribbon edge with the formation of spatially lo- 
calized breather solitons. Our molecular dynamics simu- 
lations have revealed the existence of two types of such 
breather solitons, out-of-plane vibrational breathers with 
the frequency 103 cm -1 and in-plane breathers with the 
frequency 267 cm -1 . 
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